Visualizing Atomic-Scale Negative Differential Resistance in Bilayer Graphene 
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We investigate the atomic-scale tunneling characteristics of bilayer graphene on silicon carbide using the 
scanning tunneling microscopy. The high-resolution tunneling spectroscopy reveals an unexpected negative 
differential resistance (NDR) at the Dirac energy, which spatially varies within the single unit cell of bilayer 
graphene. The origin of NDR is explained by two near-gap van Hove singularities emerging in the electronic 
spectrum of bilayer graphene under a transverse electric field, which are strongly localized on two sublattices in 
different layers. Furthermore, defects near the tunneling contact are found to strongly impact on NDR through 
the electron interference. Our result provides an atomic-level understanding of quantum tunneling in bilayer 
graphene, and constitutes a useful step towards graphene-based tunneling devices. 
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Understanding quantum tunneling at the atomic level is es- 
sential in the study of nanoscale materials and their applica- 
tions in the tunneling devices (T] |2) . One of the most in- 
teresting tunneling phenomena is negative differential resis- 
tance (NDR), characterized by the reversal of the standard 
current- voltage relationship, decreasing current with increas- 
ing voltage. NDR is the basic operating principle of Esaki and 
resonant-tunneling diodes, and has enabled various novel ap- 
plications OH]]. Motivated by its fundamental importance 
and potential applications, there have been efforts to study 
NDR in graphene 1 5 - 9 1, a prototypical two-dimensional ma- 
terial with tunable electronic properties [ |T0lfT3ll . However, 
tunneling-induced NDR has not yet been realized in graphene, 
while a recent study has proposed a new method for NDR in 
graphene that is not based on the quantum tunneling effect |9l . 

The scanning tunneling microscope (STM) is a powerful 
tool, which can not only directly probe NDR, but also provide 
key information on the mechanism. STM has been widely 
employed to study various electronic properties of graphene at 
the atomic scale, such as scattering and interference lfT4l[T5lL 
but investigation into the sub-unit cell regime has been lim- 
ited. Here, we report the first observation of NDR based 
on quantum tunneling at the vertical junction of STM over 
bilayer graphene. The applied electric field across bilayer 
graphene induces two van Hove singularities in the electronic 
spectrum, which are strongly localized on two sublattices in 
different layers. Such a localization of electronic singulari- 
ties leads to a novel atomic- scale variation of NDR, which 
is directly visualized by our high-resolution tunneling spec- 
troscopy within the single unit cell of bilayer graphene. This 
result provides the atomic-level understanding of quantum 
tunneling in bilayer graphene. 

Bilayer graphene was prepared on 6H-SiC(0001) wafers 
(N-dopant concentration of 1 x 10 18 cm -3 ) by thermal 
graphitization in a flow of argon as described in Ref. |[T6li . 
The samples were transferred through the air to the ultrahigh 



vacuum chamber (5 x 10 -11 torr), and briefly annealed up 
to 800 °C to refresh the surface. STM measurements were 
performed at 5.6 K and 78.2 K using a commercial cryogenic 
STM (Unisoku, Japan), equipped with a Nanonis controller 
(Specs, Germany). d//dV spectra were obtained by the stan- 
dard lock-in technique with a modulation voltage of 3-18 mV 
at 500 Hz and a consistent setpoint of -0.2 V and 300 pA. The 
d//dV map was taken by a grid spectroscopy, for which the 
thermal drift was carefully compensated by the atom-tracking 
method. 

The basic configuration of our experimental setup is shown 
in Fig. 1(a). It involves the tunnel junction of an STM tip 
and bilayer graphene separated by a vacuum gap of about 1 
nm. Using this technique we can measure I-V curves over a 
specific atomic site in the surface. The bilayer graphene is n- 
doped with about 0.3 eV above the charge-neutrality point 
[the Dirac energy (£d)L The interface of bilayer graphene 
with SiC(0001) produces a built-in electric field, rendering 
the charge densities in two stacked graphene layers inequiv- 
alent ifTTl . This out-of-plane symmetry breaking results in a 
bandgap with the magnitude of 0.14 eV at £b QUI - This char- 
acteristic of bilayer graphene has enabled a tunable bandgap 
either by gating lfT2l or chemical doping ifTTTl . and is also im- 
portant for NDR as explained below. 

Figure 1(b) shows an STM topographic image of bilayer 
graphene, acquired with -0.2 V sample bias and constant tun- 
neling current. The atomic structure of bilayer graphene is 
shown in the inset and overlaid with part of the image for 
comparison. Bilayer graphene consists of two Bernal stacked 
graphene layers, where the A sublattice of the top layer (A T , 
black) is located on the B sublattice of the bottom layer (Bb). 
Due to the direct coupling between them, the other sublat- 
tice in the top layer (# T , red) is dominantly imaged at low- 
bias voltage, resulting in the triangular pattern (THE). The 
full honeycomb lattice could be visualized at the higher bias 
voltage, confirming the characteristic bias dependence of bi- 
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FIG. 1: (color online), (a) Experimental configuration for the lo- 
cal tunneling spectroscopy, consisting of an STM tip and bilayer 
graphene on silicon carbide, (b) STM topographic image (1.5 x 
1.5 nm 2 ) taken on bilayer graphene with sample bias of -0.2 V and 
tunneling current of 300 pA. Inset shows the atomic structure over- 
laid with part of the STM image. Balls are carbon atoms, and solid 
(dashed) lines represent chemical bonds in the top (bottom) layer, 
(c) Point-d//dV spectrum and a corresponding I-V spectrum (inset) 
taken on the Bt site. The red line overlaid is the result of model sim- 
ulations. Red arrows highlight NDR. (d) d//dV spectra normalized 
by I-V at each initial setting current marked with different colors. 
Dotted lines in (c) and (d) indicate conductance zero. 



layer graphene (20l|2Tl. In Fig. 1(b), one can also readily 
distinguish the At site from the hollow site in the top layer 
(corresponding to Ab, blue) by their relative brightness. This 
demonstrates the high spatial resolution achieved in our mea- 
surements. 

Figure 1(c) shows a d//dV curve measured on the Bt site. 
The differential conductance value goes below zero at -0.3 
V, indicative of negative conductance (or resistance). This 
can also be directly seen in I-V measurements [inset of Fig. 
1(c)] as decreasing current while increasing the voltage from 
-0.26 to -0.34 V. This observation clearly demonstrates NDR 
in bilayer graphene. The NDR behavior was reproducibly ob- 
served in the voltage range of -0.30 ~ -0.35 V, close to . In 
Fig. 1(d), we display a series of d//dV spectra normalized by 
I IV taken with different initial tunneling currents, which are 
inversely related to the tip-sample distance (tunneling barrier 
width). They show little change and consistent NDR behav- 
iors, excluding tip-induced origins of NDR such as a local 
band bending l22l . 

An important clue for the origin of NDR can be obtained 
from its spatial distribution at the atomic scale. Figure 2(a) 
compares d//dV curves on three different atomic sites [refer 
to Fig. 1(b)]. NDR at -0.3 V occurs only on the B T site, 
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FIG. 2: (color online), (a) Point-d//dV spectra taken on three differ- 
ent atomic sites of bilayer graphene, and its spatial average offset for 
clarity. Dotted lines indicate conductance zero. (b),(c) d//dV maps 
taken with a 96 x 96 grid at the bias voltage indicated by arrows in 
(a). The inset in (b) shows the simultaneously acquired topographic 
image overlaid with the atomic structure of bilayer graphene. 



while above +0. 1 V the dominant conductance is found on the 
At site. Recording d//dV curves as a function of lateral tip 
positions, we show a differential conductance map at -0.3 V 
in Fig. 2(b). The simultaneously acquired topographic im- 
age in the inset helps determine the atomic position. Interest- 
ingly, negative conductance, represented by blue color, is cen- 
tered on each Bt atom and is surrounded by a hexagonal rim 
of positive conductance. This suggests that NDR is closely 
linked to localized states at the Bt site. On the other hand, 
the same map taken at +0.17 V shows dominant positive con- 
ductance on the A T site. This site-resolved detail has not been 
observed previously most likely due to limited spatial resolu- 
tion in the tunneling spectroscopy. Indeed, the spatial average 
of d//dV curves [the black one in Fig. 2(a)] reproduces the 
typical curve with no NDR reported previously in literature 

ED. 

To understand the origin of NDR, let us recall the band 
structure of bilayer graphene near as illustrated in Fig. 
3(a). In pristine bilayer graphene, there are two pairs of en- 
ergy bands (dotted lines) — a low-energy pair from the Bt and 
Ab atoms and a high-energy pair from the At and Bb atoms. 
When a bandgap is induced by a transverse electric field, the 
low-energy bands develop a Mexican-hat-like dispersion at 
conduction and valence band edges near the bandgap I23ll24ll . 
This rather flat dispersion leads to two distinct van Hove sin- 
gularities in the density of states | 25 ] as in Fig. 3(b). Further- 
more, these two near-gap states were theoretically predicted 
to be localized at specific atomic sites in different layers, as 
shown in Figs. 3(a) and 3(b), depending on the direction of the 
electric field (26JE3- Such a strong localization of electronic 
singularities (the flat parts of energy bands) is also important 
in understanding many-body correlations at low energies and 
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FIG. 3: (color online), (a) Band structure of bilayer graphene with 
(solid lines) and without (dashed lines) the bandgap. Red and blue 
colors in a low-energy band pair represent the atomic site, at which 
corresponding states are localized (Bt: red, Ab: blue), (b) Site- 
resolved density of states in biased bilayer graphene |26 27 1. (c) 
Energy-distribution curve at the K point taken by ARPES. Data were 
taken at the Advanced Light Source, using a Scienta R4000 electron 
analyzer (VG-Scienta, Sweden) and a 95-eV photon. Energy and an- 
gular resolutions were 30 meV and 0.1°. Red and blue arrows indi- 
cate positions of two van Hove singularities, (d) Energy diagram for 
the tunnel junction of the tip and the sample. Red solid and dashed 
lines are, respectively, the sample density of states at Vo for local 
maximum current (/max) and at Vi for local minimum current (7 m in) 
after broadening by Gaussian. The gray line is the tip density of 
states, and the shaded region is the filled tip states below Ef. 



the ground state of bilayer graphene [28], but has not been 
directly shown by experiments. 

In our case, states near -0.25 and -0.39 eV are localized 
at #t in the top layer (red) and at Ab in the bottom layer 
(blue), respectively, while states near +0.1 eV are localized 
at At in the top layer (black). Since STM is highly sensitive 
to the local density of states in the surface, this explains dom- 
inant conductance at B T for -0.2 ~ +0.1 V in Fig. 2(a) and 
that at At above +0.1 V in Figs. 2(a) and 2(c). This picture 
is fully supported by the energy-distribution curve at the K 
point [Fig. 3(c)], measured by angle-resolved photoemission 
spectroscopy (ARPES). It reveals highly asymmetric intensi- 
ties for the two van Hove singularities ifTTIl . because the pho- 
toelectrons excited from the bottom layer (Ab) are attenuated 
by scattering during their passage though the top layer. This 
also confirms that the layer localization is universal over the 
wide area covered by the photon beam (~50 fim in radius) in 
ARPES. 

The presence of a layer-polarized van Hove singularity can 



qualitatively explain NDR, as in Fig. 3(d), which shows the 
relationship between the broadened sample density of states at 
#t and the tip density of states. When the sample bias equals 
Vo, where the van Hove singularity meets the £p of the tip, 
the current reaches a local maximum, / max . As it approaches 
Vi , where the van Hove singularity no longer assists tunnel- 
ing and no state exists in the bandgap, tunneling is strongly 
suppressed, and the current decreases down to 7 m in- As the 
voltage further increases out of the bandgap, the current starts 
increasing again, resulting in the N- shaped branch in the I-V 
curve [as in the inset of Fig. 1(c)]. That is, the observed NDR 
is the result of electronic spectrum reconstructions in asym- 
metric bilayer graphene. 

For the quantitive analysis, we performed I-V sim- 
ulations based on the standard formula, I(V,z) oc 
Jg_ y DOS sample (£+y)-DOS tip (£)T(£,y,z)^, where 
DOS is the density of states and T is tunneling matrix 
element |29ti3T1l . T is estimated by the WKB approximation 
as T(E,V,z) = exp(-2 v /2(0 - V/2 - E)z, where <j> is the 
work function. We used the Gaussian-broadened local DOS 
in the theoretical calculation (26J [27) at B T with a small 
contribution of that at At (to take into account the finite 
areal averaging). The tip DOS is assumed to have a typical 
localized state of tip-apex atoms near l29l . The sample 
and tip DOS's optimized to reproduce the experimental data 
are present, respectively, by red and gray curves in Fig. 
3(d). The simulated curve [the red line overlaid in Fig. 1(c)] 
successfully reproduces the NDR behavior, confirming the 
proposed mechanism. 

To further support the proposed mechanism, we have inves- 
tigated the interaction of NDR with local defects. Recent the- 
ories predict that a single point scatterer in bilayer graphene 
develops in-gap states, suppressing the sharp van Hove singu- 
larities (32l[33l. Thus, one might expect a local suppression 
of NDR near defects. Figure 4(a) shows a large scale STM 
image of an impurity (yellow blob). The d//dV maps taken 
at two different locations about 16 nm away from defects are 
present in Figs. 4(b) and 4(c). NDR is extinguished in a sub- 
set of #t atoms, which form a periodic V3 x y/3 pattern ro- 
tated by 30° with respect to the lattice constant, and these data 
show little dependence on the origin and site of defects. The 
(a/3 x v // 3)R30° modulation is the well-known signature of 
intervalley electron scattering by a defect in bilayer graphene 
[Q31Q31. Thus, we attribute two inverse patterns in Figs. 4(b) 
and 4(c) to constructive or destructive interference of defect 
scattering, as illustrated in Figs. 4(e) and 4(f). Furthermore, 
the same map taken right next to a defect [Fig. 4(d)] shows 
that NDR completely disappears and instead strong positive 
conductance emerges on a subset of B T atoms. As a result, 
positive, zero, and negative conductance are all possible [solid 
line in Fig. 4(g)], and defects can be used to switch NDR on 
and off nearby as well as at the site more than 10 nm away 
from a defect. Such a strong impact of defects in electronic lo- 
calization should be considered for the complete understand- 
ing of scattering and interference in bilayer graphene fT4l[T5ll . 
The disappearance of NDR with the progressive suppression 



4 







-•- B1 in (b) 
B2 in (b) 
-•- B3 in (d) 







-0.5 -0.4 -0.3 -0.2 -0.1 
Sample bias (V) 



FIG. 4: (color online), (a) STM image (10 x 10 nm 2 ) with -0.05 V 
and 300 pA, where a local defect (the yellow blob) is present at the 
corner. d//dV maps taken (b),(c) at two different regions about 16 
nm away from defects and (d) at a region right next to a defect. The 
atomic structure of bilayer graphene is overlaid. (e),(f) Schematic 
illustration for two possible NDR patterns due to constructive or de- 
structive interference of defect scatterings, (g) Normalized d//dV 
spectra taken on three different Bt atoms marked in (b) and (d). Solid 
and dotted lines indicate Er> and conductance zero, respectively. 



of the van Hove singularity near the bandgap [the peak at - 
0.23 V in Fig. 4(g)] makes a clear link between them. This 
will further allow the control of NDR by tuning the bandgap 
with gating or chemical doping |[TTlfT2ll . 

Our results, combined with a few existing techniques, sug- 
gest a route towards graphene-based NDR devices. The key 
to NDR is, as demonstrated here, to exploit the near-gap van 
Hove singularity in biased bilayer graphene. In our experi- 
mental configuration, the tip can be replaced by another biased 
bilayer grahene, such that a vertical tunnel junction of two 
graphene bilayers is formed. A similar graphene heterostru- 
crture has recently been realized using ultrathin boron nitride 
films as a tunneling barrier (T][2|. In such junctions, tunneling 
would depend on sublattice-registry matching. That is, if the 
sublattices in two graphene layers facing each other, where 
the van Hove singularity is localized, were to line up, maxi- 
mal direct tunneling between van Hove singularities could be 
achieved. The tunable electronic structure of bilayer graphene 
|[TT1 l34l l35l would allow further optimization for NDR and a 
large peak-to-valley ratio of tunneling current. 
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FIG. 5: (a) Large scale STM image taken with +2 V and 50 pA, where two wide domains of monolayer and bilayer graphene are clearly 
distinguished. (b),(c) STM images at the boundary of monolayer (right) and bilayer domains (left), taken with two different bias voltages 
marked at the upper left. Insets at the upper left show magnified images of the left domains. The full honeycomb pattern at high-bias voltage 
(b) transforms into the triangular pattern at low-bias voltage (c), which is the characteristic bias dependence in bilayer graphene reported 
previously [20]. This confirms the bilayer origin of domains under study. The lower right inset in (c) shows a magnified low-bias (marked) 
image of the right domains, showing a clear honeycomb lattice of monolayer graphene. 



